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Abstract. Since the oldest globular clusters (GCs) are early residuals from the formation of the Milky Way 
(MW), GCs were exposed to the likely evolving potential of our Galaxy. The expulsion of the residual-gas from 
the GCs embedded progenitors is sensitive to the conditions in the pre-MW gas cloud. By means of N-body 
computations it is shown that gas throw-out from initially mass-segregated GCs affect the shape of the low-mass 
stellar mass function (MF) and that its imprint might still be visible in the present-day MF (PDMF). The strength 
of the tidal-field at birth influences the degree of gas-expulsion driven low-mass-star depletion and therefore the 
PDMF probes the MW potential at the time of GC formation. It is argued that among the old GC population in 
the MW, younger GCs show stronger low-mass-star loss than older GCs. This is shown to be consistent with a 
contracting and self-gravitating cloud in which fluctuations in the pre-MW potential grow with time. An initially 
relatively smooth tidal field evolved into a grainy potential within a dynamical time-scale of the collapsing cloud 
(based on 171). 



1 Introduction 



GCs of the MW may be categorized into old and young 
halo GCs (Fig.[TJi. The latter group follows an age-metallicity 
relation and is thought to have an accretion origin |2|. The 
old group of GCs is roughly coeval, spanning a metallic- 
ity range from [M/H]= -2 to -0.5. The age dispersion of 
only * 600 Myr, comparable to the free-fall time of a ho- 
mogeneous sphere with the mass and size of the MWs dark 
matter halo (1.9 X 10 12 M Q , 170 kpc, 0), is "not in con- 
tradiction with the formation from a single proto-system" 
Q], which formed the proto-MW probably along the lines 
originally suggested by Eggen et al. 1 5 1 . Therefore the MW 
GC halo consists of an in-situ (old GCs) and an accretion 
dominated, or more probably a tidally generated [6], halo 
(younger GCs) with a radial transition between the two at 
about as 10 kpc (see also [7]). The old and coeval GCs thus 
witnessed the formation of the proto-MW. 

In order to access the ambient conditions prevalent at 
GC birth we use the early violent residual-gas expulsion 
process which we show to leave an imprint in the PDMF. 
This dynamical fingerprint depends on the strength of the 
tidal-field imposed on the GCs by the proto-MW which al- 
lows us, in combination with information about the metal- 
licity (age) of the GCs, to trace back the very early events 
that contributed to the formation of the MW. 
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2 Gas expulsion and the PDMF 

Classically, computations of star clusters start with an in- 
variant IMF and the evolution of their stellar mass-function 
(MF) is determined by the evaporation of low-mass stars 
over the tidal-boundary through two-body relaxation 18191 . 
However, since clusters are observed to be segregated by 
stellar mass already at very young ages, frequently inter- 
preted as evidence for primordial mass-segregation, it may 
be questioned whether this is the whole story. Upon gas- 
expulsion and the subsequent expansion of the cluster, stars 
in the cluster outskirts are stripped off the cluster through 
the external tidal-field. If low-mass stars form at the cluster 
periphery, the low-mass star population will be preferen- 
tially depleted, leading to a flattening of the mass-function 
before classical evolution starts. 

Marks et al. [10] showed this to be the case by making 
use of the grid of residual-gas expulsion models of Baum- 
gardt & Kroupa ifTTI . If clusters expell their gas on a short 
time-scale, the stars in the cluster, assumed to be close to 
virial equilibrium before gas throw-out starts, cannot ad- 
just quickly enough to the changing potential, hence lead- 
ing to strong cluster expansion and effective loss of stars 
from the cluster outskirts. The effect of mass-loss over the 
tidal-boundary is enhanced if the star formation efficiency 
is low, i.e. a significant fraction of the total mass remains 
in the form of gas which is to be expelled, and if the tidal- 
field is strong. 

Marks et al. ifTUl show that if the cluster is mass-segregated, 
low-mass stars will be removed at first, hence flattening the 
low-mass stellar MF. If the expansion and tidal-field (and, 
thus, mass-loss) are strong, the concentration of the clus- 
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Fig. 1. Distribution of GC relative ages as a function of metallicity (taken from [4 1). Open and filled circles denote GCs with Galactocen- 
tric distance smaller and larger than 10 kpc, respectively. Clusters group into coeval old GCs (relative age of 1 corresponds tos 12.8 Gyr) 
and the younger GCs following an age metallicity relation. The age dispersion in the old group is as 0.6 Gyr, corresponding to a free-fall 
timescale (Sec.[T). Old GCs witnessed the formation of the MW. 
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concentration c=log 10 (r,/r c )' metallicity [Fe/H] 



Fig. 2. Left: Concentration parameter, c = log 10 (r,/r c ), vs. low-mass stellar PDMF slope, a, (dN/dm oc nr") in the mass-range 0.3 - 
0.8 M Q . Weakly concentrated clusters are strongly depleted in low-mass stars and no cluster with a high concentration and a depleted 
MF is found (filled dots). This trend (black solid line) can't be understood in terms of purely secular dynamical evolution. However, 
A'-body integrations of mass-segregated clusters at the time of the emergence from their birth molecular cloud (squares with error bars) 
reasonably reproduce the observed trend within the observational limits (dashed lines). Right: Global cluster metallicity, [Fe/H], vs. 
low-mass MF slope, a. Data points are labelled with the respective NGC or Pal catalogue number of the GC. Clusters having a larger 
metallicity ([Fe/H] > -1.5) are depleted in low-mass stars (a < 1) giving support to the metallicity-dependent gas expulsion scenario 
(Sec.0. 



ter decreases in parallel with the depletion of the low-mass 
star population in a mass-segregated cluster. But if expan- 
sion and the tidal-field (and mass-loss) are weak, the con- 
centration remains close to the initial value and the MF is 
only weakly depleted of low-mass stars. 

This has indeed been observed by de Marchi et al. ( IfTSll . 
Fig. [2]). That gas-expulsion is indeed at work in initiating 
this trend is strengthened by Marks & Kroupa [1], who 
show higher-metallicity GCs to have flatter present-day mass- 
functions (PDMFs). Although the observed trend had been 
recognized earlier 11 1 3 1 1 411 it lacked an explanation until re- 
cently. In particular, this trend is difficult to reconcile with 
standard dynamical evolution scenarios with a universal 
IMF as it is unclear how dynamics could possibly know 
about the metal content of a cluster. Marks & Kroupa |Q] 
interpret the PDMF trend with metallicity as evidence for 



metal-dependent cluster winds in initially mass-segregated 
clusters: As radiation couples better to metal-rich gas, the 
expulsion of the left-over gas from star formation is quicker 
and perhaps more efficient, leading to strong expansion of 
the cluster as a whole and thereby loosing preferentially 
low-mass stars over the tidal boundary. 



The de Marchi trend [12] was also unexpected from a 
classical point of view but is a natural outcome of gas ex- 
pulsion. Although it has been shown that the depletion of 
low-mass stars occurs quicker if clusters start mass-segregated 
initially 0151 . such classical scenarios without gas lack an 
understanding of the metallicity trend. 



Assembling the puzzle of the Milky Way 




Fig. 4. Contraction scenario as described in Sec. 12, II Left-hand picture: In a collapsing cloud the first clusters (the stars in the pic- 
ture) form all over the cloud in a smooth potential. The conditions are similar for all of them leading to comparable MF slopes and 
concentrations after gas expulsion. These clusters enrich their immediate environment with metals. Right-hand picture: As the collapse 
proceeds more clusters form and the potential becomes grainy enhancing the mean tidal-field strength. Star clusters, e.g. the right filled 
cluster, forming near other clusters or dense clouds experience extreme conditions. They are more strongly enriched in metals from the 
surrounding objects than other clusters, e.g. the left filled cluster which is located in a more isolated place of the grainy potential, and 
they experience stronger tidal-fields than the first-forming clusters. This leads to strong differences between the formation sites and, in 
turn, to variations in the PDMF slope. 



2.1 Assembly of the old halo of GCs 

From the models it is found that as the concentration de- 
creases, the tidal-field strength increases, i.e. present-day 
low-concentration clusters have formed in stronger tidal- 
fields (Fig. [3}. Furthermore the GCs of the old halo appear 
to be co-eval to a good approximation (Sec. [T). Despite 
their formation at the same time the clusters show a large 
spread in metallicity (-2.4 < [Fe/H] < -1.0 for the old 
halo GCs in IfTZl ). If the metal content of the more metal- 
rich old population GCs comes from the products of stel- 
lar evolution which have been produced by earlier forming 
massive stars in other clusters, the inter-cluster medium 
(ICM) will be successively enriched with time. This ma- 
terial will then be recycled into later forming clusters. In 
this case a higher metallicity as observed for the low-cr 
GCs (Fig. |2]i suggests that these clusters formed from an 
ICM enriched in metals and they should thus be somewhat 
younger than the relatively more metal-poor clusters. The 
age difference can then be at most a few hundred Myr only, 
corresponding to the free-fall time-scale (Sec.[TJ. 

These ideas lead us to a picture of the formation of the 
old inner and co-eval Galactic GC system, which might 
have formed during the contraction of the pre-MW gas 
cloud out of which finally, after an additional long period 
of ongoing accretion and merging of formerly extragalac- 
tic systems |2|, the MW has emerged (possibly along the 
lines originally proposed by Eggen et al. [5 1). These ideas 
are summarised in Fig. [4] 

In this frame, clusters which are depleted of low-mass 
stars, i.e. those that reside in stronger tidal-fields were born 
at a later stage of galaxy formation. The environment in 



which the clusters formed must then have changed dras- 
tically within a short time. The initial conditions can be 
understood if the overall potential was rather smooth in 
the beginning. The first GCs formed all over the cloud 
experiencing similar smooth tidal-fields more or less in- 
dependent of Galactocentric distance and explaining the 
comparable strengths of the tidal-fields. These first clus- 
ters enriched the ICM of their local environment with met- 
als, from which the somewhat younger clusters of the old 
population were born. 

The pre-galaxy gas cloud contracted due to self-gravitation 
during this process. The cloud may have become clumpy 
and substructures emerged. Fragmentation of the cloud into 
massive star cluster forming regions made the potential 
grainy, thereby explaining the different and on average stronger 
tidal-field strengths among the younger clusters. Clusters 
forming next to other massive, dense objects are forced to 
expel their gas under more extreme conditions (a stronger 
tidal-field) than expected from their present-day Galacto- 
centric distance. 

Gas expulsion from early forming clusters may trig- 
ger the formation of new clusters (via gas compression) 
in their immediate surrounding perhaps leading to clus- 
ter complexes, i.e. clusters of star clusters. It is imaginable 
that similar processes were at work in the GC forming gas 
cloud, increasing the average local tidal-field strengths. 

These effects lead to depleted MFs and lower concen- 
trations after gas expulsion in the younger of the old halo 
GCs. In this sense the history of events that lead to the in- 
ner GC system may involve local rapid (on time-scales of 
hundreds of Myr) re-arrangements of the interstellar mat- 
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IMF slope 
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Fig. 3. As in Fig. [2] (left) but with a modified concentration pa- 
rameter (see [ 1 1) on the abscissa. The theoretical data are coded 
for the initial tidal-field strength, in the simulations of 1 1 1 1 with 
different symbols. Open symbols and crosses depict the A'-body 
results, the filled circles labelled with their respective catalogue 
number are the data from the de Marchi GC sample. Quadrangles 
and triangles show areas with similar tidal-field strengths. All ini- 
tial models started from the dashed line, which is the slope of the 
canonical IMF (16) in the mass-range 0.3 - 0.8 M (a a 1.72), 
in which the MF slope of the observed GCs was measured. The 
tidal-field strength increases with decreasing concentration and 
increasing low-mass star depletion. Thus, clusters which expe- 
rienced stronger tidal-fields are more metal-rich (i.e. younger, 
Fig-El right panel). 

ter superposed on the overall contraction or collapse to the 
final population II spheroid. 
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